Orbital ordering in manganites in the band approach. 
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We consider the orbital ordering in LaMnOs and similar systems, proceeding from the band 
picture. We show that for the realistic magnetic structure of A-type there exists a complete nesting 
betweeen two e 9 -bands. As a result there occurs an instability towards an excitonic insulator-like 
state - an electron-hole pairing with the wave vector Q — (n, n), which opens a gap in the spectrum 
and makes the system insulating. In the resulting state there appeasr an orbital ordering - orbital 
density wave (ODW), the type of which coincides with those existing in LaMnOs. 



INTRODUCTION 

Manganite oxides Ri^M^MnOs (R=La, Pr, Nd,... 
and M = Ca, Sr) attract great interest both from ex- 
perimental and theoretical points of view. The phase 
diagram of these compounds is very rich, and in large 
part of it there exist, besides magnetic, also an orbital 
ordering due to the Jahn- Teller nature of Mn +3 ions with 
configuration t\ g e g . Usually this ordering, especially at 
the end of the series RMn03 (x=0) is described in the 
picture of localized <i-electrons p). However it was re- 
cently shown [J Q that the orbital structure of doped 
materials can be rather successfully explained from the 
opposite, band structure point of view - at least for the 
e 9 -electrons. Experimentally manganites are not very 
far from the metallic state. Although being insulating, 
undoped LaMnOs has relatively small energy gap, by 
optical data of the order of 0.4 eV 0- Moreover the 
conductivity of RMnC>3 considerably increases above the 
temperature of the orbital ordering 7] . Thus the question 
whether one can explain the main properties of undoped 
RM11O3 materials, in particular their insulating charac- 
ter and the orbital ordering, proceeding from the band 
structure point of view, has a great interest. 

The first attempt in this direction was undertaken by 
Gor'kov and Kresin Q who wrote down the spectrum of 
e g - electrons on a cubic lattice in a general form based 
only on symmetry arguments. However they did not take 
into account particular relations between hopping matrix 
elements for the e g -electrons, and did not discuss the 
specific features of the Fermi-surface, such as nesting, 
which play crucial role in our treatment 

In the present paper, following the works [MSB' 
we consider the ground state of LaMnC>3 from the band 
structure point of view. We show that in the framework 
of the double-exchange picture there exist perfect nest- 
ing between the "electron" and "hole" pockets of two e g - 
bands for the realistic magnetic structure of A-type (fer- 
romagnetic planes stacked antiferromagnetically) , and 
also for the hypothetical ferromagnetic state. As a result 



of the nesting the system is unstable against electron-hole 
pairing of the excitonic insulator type. Correspondingly 
the quasiparticle spectrum develops an energy gap, and 
the system becomes insulating. The pairing occurs with 
the wavevector Q = (tt, tt) forming two sublattices which 
correspond to orbital ordering. The type of this orbital 
ordering obtained theoretically coincides with that ob- 
served in LaMnOs experimentally, although the magni- 
tude of it and of corresponding lattice distortions in this 
weak coupling scheme is smaller than that existing in re- 
ality. Thus we show that the type of the ground state, its 
insulating character and the symmetry of orbital order- 
ing is reproduced correctly not only from the localized, 
but also from the itinerant point of view. 

We do not claim, of course, that the real LaMnC>3 is 
fully described by the weak coupling approach: in real 
systems the electron correlations arc definitely substan- 
tial. Nevertheless, the actual material may show the fea- 
tures typical for both the limiting situations, and there- 
fore we believe that the treatment carried out in this 
paper, which proceeds from the band picture, may be 
quite useful and instructive. 



THE MODEL 

In our calculations we use the tight binding approach 
for the band structure of manganites, where both dou- 
ble exchange and superexchange are incorporated. In 
the double-exchange framework with the strong on-site 
Hund's rule coupling the motion of e g -electrons is largely 
determined by the underlying magnetic structure. In 
this scheme the stability of various magnetic and or- 
bital structures is determined by the competition of the 
band energy of the e s -electrons, favoring ferromagnetism, 
and superexchange interaction J between localized (£29) 
spins, favoring antiferromagnetism. 

The effective hopping matrix elements t"- depend on 
the type of the orbitals a, (3 and relative orientation of 
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*ij|ix = etc @. Besides this, there is the usual de- 
pendence of the hopping on the angle between neigh- 
boring spins: tf 'P — t"^ cos Oij/2. The antiferromag- 
netic superexchange energy per bond is J cos % . 

Thus for this model with the infinite Hund energy we 
come to the double exchange Hamiltonian: 



(rt,Jl) (-7I,7[) 



(tt,tt) 



H — ~ E t?f C a;<ri c 0;trj 



(1) 



to which we will add later the electron-electron interac- 
tion. The localized fa g ) spins are treated classically. The 
last term in the Hamiltonian is repulsive electron-electron 
on-site interaction. 
Using 

4 = E elPRt ( a P C U-y^ + 0P c U)i) ( 2 ) 

i 

at=J2^ pRi (-^U-y^+^U)i) (3) 

i 

it is straightforward to calculate the band dispersion and 
total energy of different phases (cf. 0, ) ■ For example 
for the simplest case of A-type antiferromagnetic ordering 
one has the following bands: 



£i,2 (p) A 



COSp x + COS Py 



(4) 



± (COS 2 p x + COS 2 Py — COS Pn; COSp y ) 12 

and for the ferromagnetic case one has 



£i,2(p)A 



COS ^ + COS p y + COSp z 



(5) 



± (COS 2 p x + COS 2 p y + COS 2 Py — COSp x COS Py 



1/2 



COSp x COSp z — COS Py COS p z ) 



It is instructive to look at the form of the Fermi-surface 
for different concentration of electrons. For the A-type 
magnetic structure with increasing number of electrons n 
in e g -bands (n = 1 — x in Lai_ x Ca x Mn03) we first start 
to fill the band s\ in J2J|. It has predominantly \x 2 — y 2 )- 
character close to the T-point (center of the Brillouin 
zone). The second band £2 begins to fill from n « 0.14, 
creating the second pocket of the Fermi-surface. The 
evolution of the Fermi-surface is shown schematically in 

Fig.m 

We see that by approaching n = 1 (the case of the un- 
doped LaMnOa) the shape of the Fermi-surface pockets 
becomes more and more "flat". Although for n = 1 one 
can not speak about electrons and holes in strict sense, 
from the sequence of the filling ( Figs. 1 a-d) it is clear 
that one should treat the inner pockets of the Fermi- 
surface as electron, and the outer one - as hole-like. It is 
true at least for n slightly less than one. 
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FIG. 1: Evolution of the Fermi-surface with increasing the 
electron concentration n 
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FIG. 2: Fermi surface in the limit n — » 1. The "electron" 
pocket is shown by the diagonal hatching; "hole" - by vertical 



It is also clear that for n = 1 there exist perfect nesting 
between the bands 1 and 2: their Fermi-surfaces coincide 
with the shifting by the wavevector Q = it). It can 
be seen directly from the Eq. For n = 1 one has 



£i(p + Q) = -ea(p) 



(0) 



The spectrum for the ferromagnetic phase (3D - case) 
obeys the same relation © with the wave vector Q = 
(tt, 7t, it). Hence the full nesting between the bands 1 and 
2 take place also in the ferromagnetic case, although the 
Fermi surface is not flat anymore. 



EXCITONIC INSTABILITY AND PAIRING 

It is well known, that in the situation of nesting the 
original homogenous metallic state of the system is un- 
stable. Depending on the dominant interactions different 
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types of ordering state may exist: those with singlet or 
triplet, real or imaginary order parameters, etc. |l2( 

In order to study the consequences of nesting for mag- 
netically ordered materials we assume the magnetic order 
of localized spins of A-type. We want to note that due 
to the large Hund's coupling there is no hopping from 
one ferromagnetic plane to another and the spin indexes 
are fixed. But instead of the spin degeneracy, which is 
lifted, we have here another degree of freedom - dou- 
ble orbital degeneracy, which plays the role of spins in 
the nondegenerate case. Thus for the repulsive electron- 
electron interaction, which we assume, we would have, 
instead of the formation of a spin density wave (SDW) 
state, the corresponding formation of its orbital analog - 
orbital density wave (ODW). We will show that this state 
is indeed realized in our case. The new ground state in 
case of nesting can be easily found in the standard Green 
function technique (cf. for example 0, El). We will 
consider the case of the A-phase. The ferromagnetic case 
may be studied along the same lines. 

Due to coupling between the two bands (J2J the anoma- 
lous Green function arises. It is obvious that while nor- 
mal Green functions G 11 (p) and G 22 (p) are proportional 
to 5 P)P ', the anomalous Green function G 12 (p,p') is pro- 
portional to S p+ Qy . Further we will use notation G 12 ( P ) 
for the sake of simplicity. In general the equations for the 
Green functions can be written in the following form: 

^„-£ 2 (p)+ M )G 22 (p,c i ;„) = l + S 12 (p, Wn )G 12 (p,w„) 

(tw„ - £i(p) + fi)G 12 (p,uj n ) = E 12 (p,w„)G 22 (p,w„) 
where 



£ 12 (p,^„) = 



(7) 



r 12 (p,p',w„ - w m )G 12 (p,w m ) 



with r 12 (p,p',u; n — uo m ) being the repulsive electron- 
electron interaction. Using E\(p) = — £2(p + Q) = e(p) 
one can write: 



G 22 ( P ,^„) 



iuj n + e(p) 



^ + £ 2 (p) + |Ei 2 (p,^)| 2 

£ 12 (P,^) 
^ 2 + £ 2 ( P ) + |S]i 2 ( P ,^)| 2 



(8) 



(9) 



£ 12 ( P ,c n ) = -T]T J ^r 12 ( P , P >„- Wm ) 



£ 12 ( P ,^ m ) 



a4 + £ 2 ( P ) + |£i 2 ( P , Wm )| 2 



(10) 



Assuming T 12 ( P , P ', uj n — uj m ) being independent on fre- 
quency one gets that E 12 (p, p',w„, uj m ) is also frequency- 
independent. The quantity 



£ = ±^e(p) 2 + A 2 (p) 



(11) 



gives the spectrum of elementary excitations with the gap 
| A(p) | = |£( 12 )( P )|. The gap can be determined from the 
following equation 



a / x f dp' „ 12 , ,, A(p') 



A ^ tanh^ 



2T 



(12) 



which in the presence of the nesting always has nontrivial 
solution. For the simple on-site interaction (Hubbard 
repulsion \J n ia ^rii^ with a ^ 0) one gets in the weak- 
coupling case U <C W the usual solution with A(T — 
0) ~ £p exp(— W/U) with W being bandwidth. One can 
also show that the results do not crucially depend on 
the choice of Hubbard interaction, and that they remain 
qualitatively the same for the electron-electron repulsion 
of the screened Coulomb type with the screening length 
larger than the intersite distance. 

As we see from Eci.l[ll|l. the resulting spectrum has a 
gap and consequently this state is insulating, in agree- 
ment with the insulating nature of LaMn03. It is also 
clear that the pairing with the wave vector Q = (ir, ir) 
leads to the formation of the checkerboard-type super- 
structure in the xy-plane. 

Using the inverse transformation of J3J) one can find the 
alternation of the occupation numbers of different orbital 
states due to interband excitonic pairing (appearance of 
the anomalous average 8n^ 2 _ y2) ^ z2) = {c^ x2 _ y2)i c\ z 2 )i )) 



\ x y ) 







|x 2 -y 2 )| Z 2 > 



= (-Yd, 



(13) 
(14) 



where D = G 12 (x, x) and proportional to the energy gap. 

The obtained excitonic insulator state corresponds to 
the formation of an orbitally-ordered state in which in 
each plane there appears two orbital sublattices of the 
orbitals \<j>) = cos(0/2)|z 2 ) ± sin(</)/2)|x 2 
angle 



y 2 ) with the 



(15) 



tan(/)/2 = 

(n| z 2) - n\ x 2_ y 2))) - ^/{n\ z 2) - n^.^) 2 + 4£> 2 
2D 

where the diagonal averages n z i and n x 2_ y 2 for n = 1 



are slightly different: 



0.45, n x 2_ v 2 ps 0.55 due to 



specific form of occupied bands. For equal occupation 
n z 2 = n x 2_ y 2 — 0.5 the ordering would correspond to 
the angle (j> — tt/2, i.e. to the equal mixture of the basis 
orbitals (\tt/2 >= 1/V2(|z 2 ) ± \x 2 - y 2 ))). 

Experimentally it was found that the corresponding 
orbitals in LaMnOa are \cf> ~ ±100°); this solution in 
our analysis would be obtained for intermediate Hubbard 
interaction W/U ~ 3. We would like also to note that in 
the standard treatment of orbital ordering in the picture 
of localized electrons one obtains close to T c similar | ± 
7r/2) orbital sublattices H3)- 



4 



Turning to the ordering in the third direction, for the 
real A-type magnetic structure we see that in the ac- 
cepted approximation (no hopping and no dispersion in 
^-direction) the system "does not know" how to order. 
If we include certain weak hopping in the z-direction 
which could appear at finite temperatures and, even at 
T = 0, due to quantum corrections to the leading terms 
considered here [lfij . these hopping would lead to the 
appearance of the main instability at the wavevector 
Q = (jr, it, 7r), i.e. it would give orbital alternation also 
in the z-directions. The same type of ordering appears 
also in the ferromagnetic case. Apparently the exper- 
imentally observed ordering, which in LaMn03 is that 
with the "in phase" orbital ordering in the z-direction, 
Q = ("7r, 7T, 0), is caused by some other factors, e.g. an in- 
terplay of GdFeC>3 distortion (tilting of Mn06-octahedra) 
and Jahn- Teller deformations |l7|. 

Referring to the theoretical situation, we should note 
that in contrast to the usually considered excitonic in- 
sulators, for the A-type magnetic structure there exist 
not only an interband nesting, but for n=l also each 
pockets has nesting within itself, with its own vectors 
Qx,y = (i")0), (0, 7r) (this was also noticed by Hotta 
who, however, did not consider the interband instabil- 
ity). In principle these two instabilities may coexist and 
influence one another. 

One can show that the inclusion of both these factors 
just strengthens the excitonic instability and increases 
the resulting energy gap, without changing qualitative 
features of the solution. Besides, keeping in mind the 
sequence of band fillings shown in Figs.l a-d, we expect 
that in the realistic situation the interband pairing would 
dominate and the results presented above would be valid. 
Besides, intraband nesting exist for n = 1 for the A-type 
structure, but it is absent for the ferromagnetic case. 

CONCLUSIONS 

Thus we see that one can explain, at least qualita- 
tively, the appearance of an orbital ordering and opening 
of the gap in the electron energy spectrum of undoped 
LaMnC>3 and similar systems proceeding from the band 
picture. In this approach orbital ordering appears due to 
an excitonic instability, caused by the nesting between 
pocket of two e s -bands. The resulting state has a correct 
type and symmetry of the orbital ordering. Thus one can 
interprets the ground state of LaMnC>3 as arising due to 
this excitonic instability. 

Of course, as already mentioned in the introduction, 
we do not claim that our treatment gives the only cor- 
rect picture of LaMnOs; there are many indications that 
the electron correlations in it are relatively strong, so that 
it is usually rightly treated as a Mott insulator. However 
the truth is often "in between" , the system may show the 



features of the both the limiting pictures. In any case it 
is gratifying that one can obtain a quite reasonable de- 
scription of this system proceeding from the band point 
of view. This result, in certain sense, completes similar 
investigations carried out in 0, 0, ■ In these papers it 
was shown that the band approach works quite well for 
description of orbital structure of overdoped, half-doped 
and less-than-half-doped manganites. It is even some- 
what surprising that using the same simple approach one 
can explain even the properties of an undoped LaMn03, 
which until now was always considered as a bona fide 
Mott insulator. One may hope that this novel point of 
view may prove useful both for better understanding of 
the theoretical situation and for the interpretation of cer- 
tain experiments. 
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